Decoherence without classicality in the resonant quantum kicked rotor 
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We study the quantum kicked rotor in resonance subjected to an unitary noise defined through 
Kraus operators. We show that this type of decoherence does not, in general, lead to the classical 
diffusive behavior. We find exact analytical expressions for the density matrix and the variance in 
the primary resonances. The variance does not loose its ballistic behavior, however the coherence 
decays as a power law. The secondary resonances are treated numerically, obtaining a power-law 
decay for the variance and an exponential law decay for the coherence. 
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I. INTRODUCTION 

The development of experimental techniques has made 
possible to trap samples of atoms using resonant ex- 
changes of momentum and energy between atoms and 
laser light This progress has been accompanied by 
the development of the interdisciplinary fields of quan- 
tum computation and quantum information 0]. 

The study of open quantum systems is an outstanding 
topic of quantum mechanics. In particular the transition 
form the quantum to the classical world has intrinsic im- 
portance. On the other hand the advent of quantum 
computation makes of decoherence a central problem in 
the interaction of the quantum computer with its sur- 
roundings. 

Simple theoretical and experimental models such as 
the quantum kicked rotor (QKR) and the quantum walk 
(QW) may play an important role in this frame. Al- 
though the existent experiments have high accuracy in 
both coherent storage and manipulation of the atoms, 
the interaction with the surroundings introduces different 
degrees of decoherence influencing the unitary evolution 
of the system. 

The QKR is a milestone in the study of chaos at the 
quantum level 0] . The behavior of the QKR depends on 
whether the period of the kick is a rational or irrational 
multiple of 27r (in convenient units) [3| . For rational mul- 
tiples the behavior of the system is resonant with ballis- 
tic spreading and has no classical analog, its standard 
deviation a has the time dependence a{t) ~ t. For irra- 
tional multiples the average energy of the system grows 
in a diffusive manner for a short time and then dynami- 
cal localization take place. The quantum resonances and 
the dynamical localization of the QKR have been exper- 
imentally observed in samples of cold atoms interacting 
with a far-detuned standing wave of laser light Q and in 
particular the secondary resonances have been recently 
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observed by Kanem et al. . 

The QKR as a simple toy model allows to study the 
complexity of decoherence both analytically and numeri- 
cally, these studies have a 25 year-old history ■ On 
the other hand the first experimental observation of envi- 
ronment induced decoherence in the QKR was reported 
by Ammann et al. [lo| . 

In this line we recently investigated the QKR in reso- 
nance subjected to: a) decoherence with a Levy waiting- 
time distribution [ll|, [l^ and b) an excitation that fol- 
lows an aperiodic Fibonacci prescription [l^. In both 
cases we find that the secondary resonances have a sub- 
ballistic behavior {a{t) ~ < c < 1), while the 
principal resonances maintain the well-known ballistic 
behavior. These results are very surprising since one ex- 
pects diffusive behavior when decoherence occurs. Other 
authors also investigated the QKR subjected to noises 
with a Levy distribution [3, [l5| and almost periodic 
Fibonacci sequence [IB], showing that this decoherence 
never fully destroys the dynamical localization of the sys- 
tem, but leads to a sub-diffusion regime for a short time 
before localization appears. 

In this work we want to study the decoherence effect of 
an unitary operation described by Kraus operators (l7j 
acting on the density matrix. Our route is similar to that 
followed by Brun et. al [18] with the QW but our results 
in the QKR are very different. 

II. KICKED ROTOR 

In this section we briefly review the dynamical equa- 
tions for the QKR [3]. Its Hamiltonian is 

H ^ ^ + KcosOY^Sit-nT) (1) 

where the external kicks occur at times t — nT with n 
integer and T the kick period, / is the moment of iner- 
tia of the rotor, P the angular momentum operator, K 
the strength parameter and 9 the angular position. In 
the angular momentum representation, P\£) = £h\£), the 
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matrix element of the time-step evolution operator U is 



-ifeT/h 



J,_,(«), (2) 



where e — h^/2I, Jm is the mth order cylindrical 
Bessel function and its argument is the dimensionless 
kick strength k = K/h. The resonance condition does 
not depend on k and takes place when the frequency of 
the driving force is commensurable with the frequencies 
of the free rotor. Inspection of eq.(l2]) shows that the res- 
onant values of the scale parameter r = eT/h are the set 
of the rational multiples of 2tt, t = 2t: p/q. In what fol- 
lows we assume, that the resonance condition is satisfied, 
therefore the evolution operator depends on k, p and q. 
We call a resonance primary when p/q is an integer and 
secondary when it is not. 



III. KICKED ROTOR DYNAMICS WITH 
DECOHERENCE 

In order to generate the dynamics of the system we 
consider that the decoherence is introduced through a 
completely positive map, that is defined by a set of Kraus 
operators {An} [lB|. To preserve the trace of the quan- 
tum operation these operators satisfy 



N 



(3) 



Let us take two values of the strength parameter k: ki 
and K2. The corresponding time step operators Ui = 
U (ki) and U2 = U (^2) are used to define 



Ai = y/^ Ui, 
A2=^ U2, 



(4) 
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as a particular set of Kraus operators, where a e [0, 1] 
and [3 = 1 — a in order to satisfy Eq. ^ . Then the 
following map for the time evolution of the density matrix 
is proposed 



p{n + 1) = a Uipin)Ul + [3 U2p{n)ul, 



(6) 



where n indicates the time t — nT. When a (or /?) 
vanishes Eq.® reduces to the well known evolution of 
the usual kicked rotor in quantum resonance. In other 
cases a (or (3) may be thought as the probability per 
time-step to apply the operator Ui (or U2) to the density 
matrix. 

We shall study the previous map in the case when the 
operators Ui and U2 commute {[Ui,U2] = 0), as is the 
case in the primary resonances, in what follows we use 
the principal resonance for simplicity. In this case it is 
easy to prove, using mathematical induction, that the 
solution of the map Eq. ([6]) is 

p{n) = ui p(0) Ui^ C/^"-^"^^ (7) 



where (") = -jfj^zj^- It is important to point out that 
Eq.© is a generic solution of Eq.([6]) for any couple of 
unitary operators that commute. This means that the 
solution of the map is independent of the details of the 
model. 




FIG. 1: The coherence C(n) as a function of the dimensionless 
time n in log-log scales for Ak = 1000. The coherence was 
calculated, from top to bottom, for a = 0.1, a = 0.2, a = 0.3 
and a = 0.5. The straight stretches with slopes —0.5 show a 
power law behavior C(n) ~ 

The probability for the angular momentum value i at 
time n is P(£, n) = {l\p{n)\l). We shall calculate this 
probability for the first principal resonance. The ma- 
trix elements of Ui and U2 are expressed as {l \ Ui \j) 



Jj-i{Ki) and (;| U2 \j) 



Jj„^(k2)- Then 

using the above equation Eq. ^ with the initial condition 
p(0) = |0) (0| the probability is 



3=0 



^13' {l\U{rnmm\rn,)\l) 



(8) 

where r„j = [n — j)Ki + jK2 and (Z| J7(r„j) |0) — 
i"^ Ji{rnj)- The moments of the angular momentum are 



(£"(n)) = 



£—00 



(9) 



We want to study the time behavior of the variance 
CT^ = — (i)'^. The first moment vanishes due to the 
symmetry of the initial condition p(0). Using the prop- 
erties of the Bcssel functions, the following value for the 
variance is obtained 



{aKi + Pk2)'^ ri^ + {ki — K2y aPn . (10) 
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In the case when ki = K2 the system reduces to the 
usual kicked rotor in resonance and its variance has the 
well known ballistic behavior characteristic of this case 
When Ki ^ K2 the coherence of the system is lost, 
as it is shown below, because the probabilistic map is 
effectively working. Eq. (flO|) leads us to some interesting 
results. It shows that the ballistic behavior is maintained 
with this decoherence; but additionally there appears the 
diffusive term (ki — ^2)^ a(3n. In particular if the param- 
eters verify aKi + (3k2 — 0, the behavior of the variance 
is totally diffusive as in the classical random walk. Then 
we can conclude that this decoherence always affects the 
behavior of the variance but, in general, does not break 
its ballistic growth. 

The degree of coherence of the system can be measured 
by several means. We choose the following: 

1=0 

Substituting Eq.([7]) in the above equation, and using the 
properties of the Bessel function the equation for the co- 
herence is obtained 



j=0 i=0 



'73^a"-^/3V2(AA^,,), (12) 



where Anij — (z — j) Ak, with Ak = ki — K2- From 
Eq.lHl]) is easy to prove that C(0) = 1 and C{n) < 1 
for n > 0, but in general this equation will be difficult 
to reduce to a more simple expression. However we can 
get some additional information when An is very large. 
In this case Jq (Ak^ ) goes to zero, except when i = j 
because Jq (O) = 1. Then in this limit Eg. p^ reduces to 



j=o 



(13) 



Here we observe the interesting result that the coherence 
is independent of the strength parameters of the system 
if Ak is sufficiently large. 

We made numerical studies of the long-time behavior 
for Eq. (fT2)) as a function of the parameter a. In Fig. [T] 
the function C{n) is plotted for large Ak and different 
values of a. This figure shows a power law decay for the 
coherence C{n) ~ independently of the value of a. 

The same results were obtained using Eg. p^ . Therefore 
in Eq. p3|) we may choose a particular value of a to cal- 
culate its long-time decay. Taking a = 1/2 and using the 
sums of the binomial coefficients we obtain 



Cin) 



2n 



(14) 



For large n is possible to use the Stirling formula to ob- 
tain analytically the following expression for the coher- 
ence 



C{n) 



(15) 



confirming our numerical result. 

We also studied the coherence for several smaller values 
of Ak. We obtained that for each fixed value of Aft the 
coherence always decays as the power law C(n) ~ 
with 7 > (see Fig. [2]). Additionally we observed that 
the exponent 7 is always independent of a like in Fig. [U 
therefore 7 only depends on Ak. We observe that the 
exponent 7 grows with Ak with its values in the interval 
[0.3,0.5]. We can conclude that the qualitative behavior 
of Eq. (fT2| and Eq. (fT3|) are the same for all values of Ak. 



1 - 



0.01 




10 



100 

n 



~ I I I I I 

1000 



FIG. 2: The coherence C{n) as a function of the dimensionless 
time n for a — 0.1 in log-log scales. For large n the curve 
satisfy a power law C(n) ~ n~^. The parameters of the 
curves, from top to bottom, are (1) Ak = 0.2 and 7 — 0.36, 
(2) Ak = 0.3 and 7 = 0.37, (3) Ak = 1 and 7 = 0.4, (4) 
Ak = 10 and 7 = 0.47 



Now we inquire the incidence of decoherence on the 
secondary resonances. In this case the commutativity 
between the evolution operators Ui and U2 is lost and 
the expressions for the variance and the density matrix 
become very cumbersome. Then we study the decoher- 
ence numerically using Eq.® for several values of the 
parameters k and a. 

In Fig. [3] the standard deviation a is presented, for 
fixed values of ki and K2 and for different values of a. It 
is seen that a has power-law decay with an exponent 
c that depends on a. This parametric dependence is 
very different from that given by Eg. ljlOp in the primary 
resonances where c = 1. The values of c were adjusted 
for the last thousand values of n and we found that they 
are near c = 1. For other values of ki, K2 and a the 
exponent c varies between 0,4 and 1.2. Considering all 
the cases studied we conclude that the exponent c does 
not show a clear rule of dependence with the parameters. 

The numerical study of the coherence C(n), for the 
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FIG. 3: The standard deviation cr(n) as a function of the 
dimensionless time n. The parameters are ki = 0.1, K2 — 0.2, 
p/q = 1/3 in all cases. Dashed line a — 0.1 (c = 1.2), thick 
line a = 0.2 (c = 1.), thin line a = 0.5 (c = 0.9). 

same range of values of the parameters as for a, showed 
that its time decay is better approximated by an expo- 
nential than by a power law, i.e. C{n) ~ exp {—Sn) with 
(5 > 0. Therefore the coherence of the system in the 
secondary resonances is lost faster than in the primary 
ones. 



IV. DISCUSSION AND CONCLUSION 

Decoherence in quantum systems as QKR or QW has 
been extensively studied. Analytical and numerical re- 
sults [H, [H, [13, mi on the effect of different kinds of 
noise have shown that quantum properties are highly sen- 
sitive to random events. In particular the linear increase 
of the standard deviation a{t) t can be eventually sub- 
stituted by a diffusive behavior <7{t) ^ t^l"^ as in the 
classical random walk. 

The linear increase of the standard deviation of the 
QKR in resonance is usually accepted as a direct con- 
sequence of the quantum coherence, in other words, a 



consequence of the unitary evolution. This work shows 
explicitly that unitary decoherence does not break the 
temporal linear increase of a. 

The absence of diffusive behavior in presence of deco- 
herence has already been shown in our previous works 
[Til \vX [Tsj. There we have studied the QKR subject to 
different types of noise with a Levy waiting time distri- 
bution and we found that the system has a sub-ballistic 
wave function spreading and its standard deviation has a 
power-law tail. However in that opportunity the coher- 
ence had not been studied. 

Here we have considered a new type of decoherence in 
the QKR as an unitary map acting on the density matrix. 
We obtain an analytical expression for the density matrix 
when the Kraus operators commute. We prove that the 
decoherence affects the variance but its ballistic growth 
persists in spite of an additional linear term. Therefore 
asymptotically the linear behavior of the standard devi- 
ation is not suppressed by the noise. On the other hand 
the coherence C(ri) has a power-law decay for all values 
of the parameters. We want to underline that the den- 
sity matrix Eq.(l7]), solution of Eq.® only depends on 
the commutativity of the unitary operators Ui, U2 and it 
is independent of their detail. This allows to extend the 
use of this expression for other quantum models such as 
the QW. In previous works [ll[, [ij, [H, [ll] we have estab- 
lished a parallelism between the QKR in resonance with 
the discrete QW. Then the type of treatment presented 
in this paper could be applied to the QW. 

When the Kraus operators do not commute we have 
not usable analytical expressions, it is necessary to make 
numerical studies. We establish that: a) the standard 
deviation has no simple dependence with the parameters 
of the system, b) the standard deviation has (in the long- 
time limit) a continuous range of behaviors from diffusive 
to ballistic and c) the coherence C (n) shows a exponential 
law decay. 

We can conclude that the effect of decoherence of the 
type studied in this work does not necessarily transform 
our quantum system into a dissipative system such as a 
Markov process. In more general terms, the mere pres- 
ence of noise does not assure the passage from the quan- 
tum to the classical world. 
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